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Naively resummed perturbative approximations to the thermodynamic functions of QCD 
do not converge at phenomenologically relevant temperatures. Here we review recent results 
of a three-loop hard-thermal-loop perturbation theory calculation of the thermodynamic 
£\J , functions of a quark-gluon plasma for general N c and Nf. We show comparisons of our re- 

cent results with lattice data from both the hotQCD and Wuppertal-Budapest groups. We 
demonstrate that the three-loop hard-thermal-loop perturbation result for QCD thermody- 
namics agrees with lattice data down to temperatures T ~ 2T C . 
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§1. Introduction 



A key question in the study of the quark-gluon plasma (QGP) is whether or not 
one can use a weakly-coupled framework to calculate properties of the system at tem- 
peratures near the QGP phase transition temperature, T c ~ 170 MeV. This is phe- 
nomenologically relevant because the ultrarelativistic heavy-ion collision experiments 
at Brookhaven National Labs (RHIC), and forthcoming at CERN (LHC) generate 
initial QGP temperatures which are on the order of 2T C and 5T C , respectively. At 
these temperatures the strong coupling constant is approximately a s = g 2 /47r ~ 0.3. 
Unfortunately, this places the system in regime where the coupling is neither very 
small nor very large. Prior to experiments at RHIC, theorists expected that the 
QGP could be described in terms of weakly interacting quasiparticles; however, data 
from RHIC suggested that the state of matter created behaved more like a strongly 
coupled fluid with a small viscosity.—' This has inspired work on strongly-coupled 
formalisms based on e.g. the AdS/CFT correspondence. However, some observables 
such as jet quenching^ and elliptic flow® can also be described using perturbative 
methods and so it is difficult to decide whether the plasma is strongly or weakly 
r> | coupled based only on RHIC data. 

In the upcoming heavy-ion collisions at LHC, the energy densities and therefore 
the initial temperatures will be higher than those achieved at RHIC. One expects 
initial temperatures on the order of 5T C . An important question is then whether 
the matter generated can be described in terms of weakly interacting quasipar- 
ticles at these higher temperatures. Lattice simulations of QCD provide a clean 
testing ground for the quasiparticle picture and in this proceedings report we com- 
pare recently obtained next-to-next-to-leading order (NNLO) results for thermody- 
namic functions of QCD 4 ' with lattice dataP® and with previous results at lead- 
ing order (LO) and next-to-leading order (NLO).^ The calculations presented here 
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are based on the hard-thermal-loop perturbation theory (HTLpt) reorganization of 
finite-temperature perturbation theory. In HTLpt one expands around an ideal 
gas of massive gluonic and quark quasiparticles where screening effects and Landau 
damping are gauge-invariantly built in. Our results indicate that lattice data are 
consistent with a quasiparticle picture down to temperatures of T ~ 2T C . 

The free energy of QCD is now known up to order log a s P® Unfortunately, 
a straightforward application of perturbation theory is of no quantitative use at 
phenomenologically relevant temperatures. The problem is that the weak-coupling 
expansion oscillates wildly and shows no sign of convergence unless the temperature 
is extremely high. For example, if one compares the ^-contribution to the QCD 
free energy with three quark flavors to the ^-contribution, the former is smaller 
only if a s < 0.07, which corresponds to T ~ 10 5 GeV or T ~ 5 x 10 5 T c . There are 
several ways of reorganizing the perturbative series at finite temperature.^ 1 These 
reorganizations are based on a quasiparticle picture where one is perturbing about 
an ideal gas of massive quasiparticles, rather than that of an ideal gas of massless 
quarks and gluons. In scalar ^-theory the basic idea is to add and subtract a thermal 
mass term from the bare Lagrangian and to include the added piece in the free part 
of the Lagrangian. The subtracted piece is then treated as an interaction on the 
same footing as the quartic termPJ In gauge theories, however, simply adding and 
subtracting a local mass term, violates gauge invarianceP^ 1 Instead, one adds and 
subtracts an HTL improvement term, which dresses the propagators and vertices 
self-consistently so that the reorganization is manifestly gauge invariant.^ 1 



§2. Hard-thermal-loop perturbation theory 



The Lagrangian density for an SU(iV c ) Yang-Mills theory with Nf fermions in 
Minkowski space is 

£qcd = -^Tr [G^CT] + ty-fD^ + £ gf + C gh + AC QC d , (2-1) 

where the field strength is G^ u = d^A u — d v — ig s [A^,A u ] and the covariant 
derivative is = — ig s A^. ACqcd contains the counterterms necessary to 
cancel the ultraviolet divergences. The ghost term £ g h depends on the gauge-fixing 
term C g f. In this paper we choose the class of covariant gauges where the gauge- 
fixing term is £ g f = — |Tr (d^A^) 2 . With the standard normalization, we have 

c A = N c , d A = N 2 -l,s F = N f /2, d F = N c Nf, and s 2F = {N 2 C - l)N f /AN c . 

Hard-thermal-loop perturbation theory is a reorganization of the perturbation 
series for thermal QCD. The Lagrangian density is written as 

C = (C Q cd + C HTL ) +AC H tl, (2-2) 

where Z\£htl contains the additional counterterms necessary to cancel the ultravi- 
olet divergences introduced by HTLpt. The HTL improvement term is 

1 '- " ,2ubi iJ iJ \ nu, \ . /-, e\- 27 W y 



Attl = -2^-^ m D^[ G ^\^hyj G %\ + (l-*)*™?^ 1 V>, (2-3) 
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where y^ = (l,y) is a light-like four-vector, and (. . .) y represents the average over 
the directions of y. The free parameters mo and m q are identified with the Debye 
screening mass and the fermion thermal mass, respectively. The parameter 5 is 
a formal expansion parameter and bookkeeping device: HTLpt is defined as an 
expansion in powers of 5 around 5 = 0. This expansion generates systematically 
dressed propagators and vertices. It also automatically generates new higher-order 
terms that ensure that there is no overcounting of Feynman diagrams. We emphasize 
that HTLpt is gauge invariant at each order in 5. 

The HTL perturbative expansion generates new divergences. We use MS di- 
mensional regularization with a renormalization scale fi to regularize infrared and 
ultraviolet divergences. There is no rigorous proof that HTLpt is renormalizable, so 
the general structure of the counterterms is unknown. However, one can show that 
through NNLO, HTLpt can be renormalized using only local counterterms for the 
vacuum, the Debye and fermion masses, and the coupling constant. The counterterm 
for a s coincides with the perturbative value giving rise to the standard one-loop run- 
ning. We do not list the counterterms here, but present the full results elsewhere.^ 
If the expansion in 5 could be carried out to all orders, the final result would be 
independent of the HTL parameters mr> and m q . However, at any finite order in 
5, the results depend on mrj and m q . A prescription is then required to determine 
these parameters. We will discuss the prescription we use below. 



§3. Thermodynamic potential 



In this section, we present the final results for the thermodynamic potential 
Q at orders 5° (LO), 6 (NLO), and 5 2 (NNLO). The LO and NLO results were 
first obtained previously^ and they are listed here for completeness. At LO, the 
thermodynamic potential was calculated exactly, while at NLO and NNLO the re- 
sulting expressions for the diagrams are too complicated. To make the calculations 
tractable, the thermodynamic potential is therefore evaluated approximately by ex- 
panding them in powers of mu/T and m q /T which assumes that these ratios are 
0{g s ). This implies that the thermodynamic potential is evaluated in a double ex- 
pansion in g s , uid/T, and m q /T, and we have kept terms that contribute naively 
through order g^. Due to the magnetic mass problem,^ HTLpt suffers from the 
same infrared divergences as ordinary perturbation theory and g^ is the highest order 
computable using only perturbative methods. 

The complete expression for the leading order thermodynamic potential is given 

bj0 

^lo -. . 7 d F 15. 2 d F . 2 , gn -3 , 45 /, A 7 , , 71-2 



60^(vr 2 -6)rf4 , (3-1) 



where Jideal = ~{^c ~ l)7r 2 T 4 /45 is the free energy of an ideal gas of noninteracting 
gluons and is the Euler-Mascheroni constant. Moreover, we have introduced the 
dimensionless parameters fi = fi/2irT, rhrj = mr)/2TrT, and m q = m q /2irT. 
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The NLO thermodynamic potential reads 7 -* 
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Finally, our new result for the NNLO thermodynamic potential for QCD is 
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where £(z) is the Riemann zeta-function. 
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Fig. 1. Comparison of LO, NLO, and NNLO predictions for the scaled pressure for Nf = 2 + 1 (left 
panel) and Nf =2 + 1 + 1 (right panel) with lattice data from Cheng et al.^ and Borsanyi et 
alP We use N c = 3, three-loop running for a s , /i = 2nT, and Aj^ = 344 MeV. Shaded band 
shows the result of varying the renormalization scale p by a factor of two around /i = 2ixT for 
the NNLO result. See main text for details. 



As pointed out earlier, the HTL mass parameters are completely arbitrary and 
we need a prescription for them in order to complete a calculation. The variational 
mass prescription unfortunately gives rise to a complex Debye mass and m q = at 
NNLO. One strategy is therefore to throw away the imaginary part of the thermody- 
namic potential to obtain thermodynamic functions that are real valued QSJE) Here 
we use another strategy explored in Refs. [T6]) . ll7p that is inspired by dimensional 
reduction: We equate the Debye mass with the mass parameter of three-dimensional 
electric QCD (EQCD),3 i.e. mp = the- In Ref. 9), it was calculated to NLO giving 



2 4vra s 2 J c\a s ( 5 11 11 p\ c A s F a s 
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This mass can be interpreted as the contribution to the Debye mass from the hard 
scale T and is well defined and gauge invariant order-by-order in perturbation theory. 
However, beyond NLO, it will also depend on factorization scale that separates the 
hard scale and the soft scale gT. For the quark mass, here we choose m q = 0. 

The final NNLO results are very insensitive to whether one chooses a perturba- 
tive mass prescription for m q or m q = 0; however, convergence is improved with the 
choice m q = 0. A detailed presentation of the full calculation of the NNLO thermo- 
dynamic potential and the dependence of our final results on the mass prescriptions 
for rri£) and m q is forthcoming in a longer paper P3J 

§4. Results 



In Fig. [H we show the normalized pressure for N c = 3 and Nf = 2 + 1 (left 
panel), and N c = 3 and Nf = 2 + 1 + 1 (right panel) as a function of T. The results 
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at LO, NLO, and NNLO use the BN mass given by Eq. (j3-4[) as well as m q = 0. 
For the strong coupling constant a s , we used three- loop running^ with A^s = 344 
MeV which for Nf = 3 gives a s (5 GeV) = 0.2034PJ The central line is evaluated 
with the renormalization scale = 2irT which is the value one expects from effective 
field theory calculations®'^ 1 and the band represents a variation of fi by a factor of 
two around this scale. 

The lattice data from the Wuppertal-Budapest collaboration uses the stout ac- 
tion and have been continuum extrapolated by averaging the trace anomaly measured 
using their two smallest lattice spacings corresponding to iV T = 8 and N T = 1OP0 
Using standard lattice techniques, the continuum-extrapolated pressure is computed 
from an integral of the trace anomaly. The lattice data from the hotQCD collabora- 
tion are their N T = 8 results using both the asqtad and p4 actionsP* The hotQCD 
results have not been continuum extrapolated and the error bars correspond to only 
statistical errors and do not factor in the systematic error associated with the cal- 
culation which, for the pressure, is estimated by the hotQCD collaboration to be 
between 5 and 10%. We note that there are hotQCD results for physical light quark 
masses; 21 ' however, these are available only for temperatures below 260 MeV and 
the results are very close to the results shown in the figures so we do not include 
them here. 

As can be seen from Fig. [T] the successive HTLpt approximations represent an 
improvement over the successive approximations coming from a naive weak-coupling 
expansion; however, as in the pure-glue casej 1 -^ the NNLO result represents a sig- 
nificant correction to the LO and NLO results. That being said the NNLO HTLpt 
result agrees quite well with the available lattice data down to temperatures on the 
order of 2T C ~ 340 MeV for both Nf = 3 (Fig. □ left) and N f = 4 (Fig. □ right). 
Below these temperatures the successive approximations give large corrections with 
the correction from NLO to NNLO reaching 100% near T c . 

In Fig. [21 we show the NNLO approximation to the trace anomaly (interaction 
measure) normalized to T 4 as a function of T for 7V C = 3 and Nf = 3 (left panel) 
and for N c = 3 and Nf = 4 (right panel). In the left panel we show data from 
both the Wuppertal-Budapest collaboration and the hotQCD collaboration taken 
from the same data sets displayed in Fig. [1] and described previously. In the case 
of the hotQCD results we note that the results for the trace anomaly using the 
p4 action show large lattice size effects at all temperatures shown and the asqtad 
results for the trace anomaly show large lattice size effects for T > 200 MeV. In the 
right panel we display a parameterization (solid blue curve) of the trace anomaly for 
Nf = 4 published by the Wuppertal-Budapest collaboration® since the individual 
data points were not published. In both the left and right panels we see very good 
agreement with the available lattice data down to temperatures on the order of 
T ~ 2T C . 



*' We note that the Wuppertal-Budapest group has published a few data points for the trace 
anomaly with N T = 12 and within statistical error bars these are consistent with the published 
continuum extrapolated results. 
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Fig. 2. Comparison of NNLO predictions for the scaled trace anomaly with Nf = 2 + 1 (left panel) 
and Nf = 2+1 + 1 fermions (right panel) lattice data from Cheng et al.^ and Borsanyi et alP^ 
We use N c = 3, three-loop running for a s , fi = 2nT, and Aj^q = 344 MeV. Shaded band shows 
the result of varying the renormalization scale ^ by a factor of two around fi = 27rT. See main 
text for details. 



§5. Summary and outlook 

We have reviewed recent results for the LO, NLO, and NNLO thermodynamic 
functions for SU(-/V C ) Yang-Mills theory with Nf fermions using HTLpt. We com- 
pared our predictions with lattice data for N c = 3 and Nf 6 {3, 4} and found that 
HTLpt is consistent with available lattice data down to T ~ 2 T c for the pressure and 
the trace anomaly. This is in line with expectations since one is expanding about 
the trivial vacuum = and therefore neglects the approximate center symme- 
try Z(N C ). Close to the deconfinement transition, it is essential to incorporate this 
symmetry P3 

Comparing our results with the NNLO results of pure Yang-Mills,^- we find 
that including the quarks gives much better agreement with lattice data. This is not 
unexpected since fermions are "perturbative" in the sense that they decouple in the 
dimensional-reduction step of effective field theory. 

As was the case with pure Yang- Mills we found that the variational solution for 
the Debye mass mc is complex and we therefore chose instead to use the perturbative 
mass parameter from EQCD together with m q = 0. Whether the complexity of the 
variational Debye mass is due to the additional expansion in mu/T and m q /T is 
impossible to decide at this stage. We also found that there was a large correction 
going from NLO to NNLO. Unfortunately, due to the magnetic mass problem it is 
impossible to go to N 3 LO to see whether the problem persists without supplementing 
our calculation with input from three-dimensional lattice calculations. 

In closing, we emphasize that HTLpt provides a gauge invariant reorganization 
of perturbation theory for calculating static and dynamic quantities in thermal field 
theory. Given the good agreement with lattice data for thermodynamics, it would be 
interesting to apply HTLpt to the calculation of real-time quantities at temperatures 
that are relevant for LHC. 
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